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We report families of discrete optical solitons in frequency space, or spectral-discrete solitons existing in a 
dispersive Raman medium, where individual side-bands are coupled by coherence. The associated time-domain 
patterns correspond to either trains of ultrashort pulses, or weakly modulated waves. We describe the physics 
behind the spectral localization and study soliton bifurcations, stability and dynamics. © 2008 Optical 
Society of America 

OCIS codes: (190.5650) Raman effect, (060.5530) Pulse propagation and solitons 



Discrete solitons in systems of coupled waveguides 
have been a subject of intense recent interest, due to 
the richness of optical effects and the potential appli- 
cations associated with them. 1,2 A seemingly different, 
and also very active subfield of optics is the generation of 
ultrashort pulses via the excitation of multiple harmon- 
ics in Raman-active and other materials. 3 In this work, 
we reveal and explore a surprising link between discrete 
optical solitons and short pulse generation in a Raman 
medium. 

It can be seen from numerical results 4, 5 that the spec- 
tral harmonics involved in short-pulse generation by Ra- 
man coherence, can either be spread over the entire fre- 
quency space covered by the model equations, or be ac- 
cidentally localized around some central frequency. The 
effect of spectral localization has not attracted much at- 
tention so far and it is often assumed that cascaded gen- 
eration of Raman side-bands is practically limited by 
only the frequency dependence of Raman gain or ma- 
terial losses. Below, we neglect the above two effects 
and demonstrate the existence of spectrally localized 
structures supported by a balance between the normal 
group velocity dispersion (GVD) of a Raman medium 
(which, quite surprisingly, plays the role of a holding 
potential) and the nonlinear coupling between the Ra- 
man side-bands (playing the role of the effective discrete 
diffraction happening in the frequency space). These lo- 
calized structures bear some features of discrete soli- 
tons in waveguide arrays, 1,2 but at the same time have 
many unique properties. It appears natural to call them 
spectral-discrete solitons (SDSs). Similar ideas of explor- 
ing localization in Fourier space have been recently used 
in a different physical context for an interpretation of 
the Fermi-Pasta-Ulam paradox. 6 

We start with the well-established model 3-5 describing 
the evolution of 2N + 1 side bands in a Raman medium 
coherently excited away from resonance. In the dimcn- 
sionlcss form, the corresponding equations are 



id z E n = (3 n E n + q*E n -i + qE, 



n+l, 



(1) 



Here n = — N, . . . , 0, . . . , N, z is the dimensionless prop- 
agation coordinate and E n are the amplitudes of the 
side-bands with frequencies u) n = loq +n(u>R + jn). ljr is 



the Raman transition frequency and j[i is the detuning 
of the frequency difference of the neighboring harmonics 
from ujr. It has been convenient for us to split this de- 
tuning into its sign j = ±1 and its absolute value [i > 0. 
P n is the propagation constant at the frequency u n . The 
frequency dependence of the propagation constant (3(ui) 
can be fitted with a polynomial in lu, which in the dis- 
cretized frequency space is 
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where p\ is associated with the group velocity at n = 0, 
P2 with the group velocity dispersion, and all the higher 
order coefficients correspond to the higher order disper- 
sions. 

q is the Raman coherence responsible for coupling be- 
tween the side-bands and ensuring their cascaded excita- 
tion. Any possible frequency dependence of the coupling 
is neglected for simplicity, q is expressed as a product of 
the amplitudes ai,2 (|ai| 2 + \ a 2\ 2 = 1) of the states sep- 
arated by lur: q = 2a\a\. The normalized Hamiltonian 
for the two-photon Raman process (with the decay and 
Stark terms omitted) is given by 3-5 
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S = ^ E nE* 



(3) 



If S varies slowly in comparison to the separation of 
the H eigenvalues, then from idta = Ha, where a = 
(ox, d2) T , an approximate expression for q can be found: 
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where j = +1/ — 1 corresponds to the in-phase/anti- 
phase dressed state, respectively. 3, 5 

If q is real, the correspondence between Eqs. Q and 
the equations describing propagation in an array of cou- 
pled linear waveguides 1,2 is obvious: q plays the role of 
the coupling between the waveguides. However, in our 
case the coupling itself is a nonlinear function of the 
field amplitudes, see Eqs. (@J, ©■ The combination of 
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equations Q and 10} is invariant with respect to the 
symmetry transformation 

E n ^E n e 1 ^, (5) 

where ip is an arbitrary constant. Thus the phase dif- 
ference between the neighboring side-bands is an arbi- 
trary number, which does not change the properties of 
any solution found. For q real, its sign, s = sign(q), is 
controlled by the relative phases of the Raman harmon- 
ics. The in-phase neighboring harmonics correspond to 
ijj = and s — 1, while the out-of-phasc harmonics im- 
ply that "0 = 7r and s = —1. However, the overall sign of 
q*E n -i + qE n+ i is not changed under the transforma- 
tion J5j) with any ip, including tp = tt. This is because the 
phase of q, see Eq. Q, compensates the phases of E n ±\. 
It contrasts sharply with the discrete waveguide array 
model, 2 where Eq. (J5J) is not a symmetry, and therefore, 
taking ip = tt reverses the discrete diffraction sign, which 
strongly effects propagation in the array. 1,2 

It is also useful to consider the continuous limit with 
E n -> E(u), fin -> and E n±1 -> E(w)±6wd u E + 
\(5us) 2 ■ d^E + 0(Slu 3 ), where 8u = (ujr + jn)/uj - In 
this limit, Eqs. QJ, Q transform into the Schrodingcr 
equation for E(w): 

ia,E = t3(u)E + j\q\^ + 2j\q\E, (6) 

the /3(w)-term forms a potential in frequency space, while 
q is the coefficient in front of the effective diffraction term 
controlling the spreading of a wavepacket in frequency 
space. V is the total power, which is a conserved quantity 
(d z V = 0) and so in the continuous limit, q itself is 
a constant. The third term in © is simply a constant 
shift of the propagation constant. We note that j = —1 
corresponds to the sign in front of the second derivative 
being that which is usual in the quantum mechanical 
context. Therefore, if (3{ui) has a minimum at some uj 
it indicates the presence of a potential well in frequency 
space. 

The simplest parabolic potential well is formed by 
Pi > and Pk>2 = 0, corresponding to the normal GVD 
typical for Raman gases. 3 It is clear that by taking ac- 
count of GVD only, i.e. Pk>2 = 0, we reduce Eq. (jSJ 
to the harmonic oscillator problem, which has a known 
set of localized eigenmodes. Thus Eqs. (Q, (0} are also 
expected to have a set of localized solutions. In a way, 
this is analogous to the approach used for fiber solitons. 
Their existence was first discovered theoretically in the 
approximation of zero higher-order dispersions, and later 
their relative robustness to the realistic dispersion pro- 
files was verified. 2 To find stationary solutions to the 
nonlinear system QJ, |J3} we make the ansatz 

E n (z) = B n exp(iKz + ianz) (8) 

and numerically solve the resulting algebraic equations 
for B n using a Newton method. K and a are the two 




Fig. 1. Single-peak (ground state) and two-peak SDSs 
calculated for j = — 1 . k varies from —0.004 to —0.93 in 
(a) and from 0.009 to —0.79 in (b), which gives \q\ from 
0.005 to 0.5. p 2 = 0.01, fi=l,a = 0, p k >2 = 0. 




K 

Fig. 2. Solid lines show power of the first five SDSs vs 
\q\ (a) and vs k (b) for j = —1. Dashed lines in (b) 
show j = +1 SDSs. The other parameters as in Fig. 1. 
Black/gray lines mark stable/unstable solutions. 

parameters enabled by phase symmetries of our system. 
We notice that substitutions K — > k — po and a — > a — 
Pi eliminate parameters po,i- In other words, without 
restriction of generality one can assume po,i = in Eq. 

Using the above approach we have been able to find 
families of spectrally localized solutions, which we called 
spectral-discrete solitons (SDSs). Examples of SDSs cor- 
responding to the first and second modes of the harmonic 
oscillator are shown in Figs. 1(a), (b) as functions of 
the coherence \q\. Doing numerical continuation of SDSs 
along k we computed q at each step and used \q\ as a pa- 
rameter in Figs. 1,2. |g| tends to its maximum, \q\ = 1/2, 
when \S\ 3> 1. Notably, the derivative of \q\ in \S\ tends 
to zero under the same conditions, i.e. |g| becomes field 
independent. Thus, perhaps counter- intuitively, in the 
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Fig. 3. Propagation of the unstable two-peak SDS with 
noise added (a) and with its power increased by a regular 
correction (b). K = 0.005 and \q\ « 0.013. The other 
parameters as in Fig. 1. 

limit of the maximum Raman coherence Eqs. (|TJ become 
quasi-linear. Oppositely, the nonlincarity of our system 
is at its maximum for |5| <C 1. The dependencies of the 
power P = ^2 n \E n \ 2 of the first five SDSs (labeled as 
T1-T5 in all Figs.) as functions of both \q\ and K are 
shown in Figs. 2(a,b). Black lines in Fig. 2 correspond 
to the stable SDSs and gray ones to the unstable ones. 
The boundaries of existence in k of some of the SDSs 
(see vertical lines in Fig. 2(b)) are given by the linear 
spectrum of Eqs. 0J. 

\q\ tends to 1/2, when k is negative and its absolute 
value increases. For these parameter values, the shape of 
SDSs is very close to the modes of the linear harmonic 
oscillator. Since nonlincarity is diminished under these 
conditions, it is not surprising that SDSs are stable for 
k < 0. Diffraction in the frequency space increases for 
large and so is the number of excited Raman har- 
monics. Thus in the time domain, the quasi-linear SDSs 
represent a periodic sequence of ultra-short pulses, see 
the full line in the inset to Fig. 2(a), where the intensities 
of the total field / = E n (z) exp(ioj n t)\ are plotted. 
Here, t is dimensionless time. At intermediate and low 
values of Raman coherence, the diffraction in frequency 
space is relatively small and nonlinear interaction starts 
to play an important role. Both of the above factors 
contribute to stronger spectral localization of SDSs, see 
Fig. 1. The corresponding temporal profiles transform 
into weakly modulated wavetrains, see the dashed curve 
in the inset to Fig. 2(a). Dynamical instabilities of the 
higher order SDSs are readily found in this regime, see 
Figs. 2(a,b). However, the ground state solution remains 
stable within the entire range of its existence. An exam- 
ple of evolution of spectral components of the unstable 
two-peak SDS perturbed by noise is shown in Fig. 3(a). 
One can see that localization due to potential (3{uj) re- 
mains unaffected by the instability development. Shape 
deformation of the unstable SDS increasing its power 
leads to formation of spectral breathers, see Fig. 3(b). 
Note that the dependencies of P on k in Fig. 2(a) indi- 
cate that transition from the stable to unstable regimes 
may be driven by the Vakhitov-Kolokolov type of insta- 
bility. 2 

The above discussion focuses on the SDSs branches 
found for negative Raman detuning j = — 1 (T-type 



Fig. 4. (a) An example of SDS for j = +1. k = 0.01, \q\ w 
0.00035 (black lines and circles), K = 0.02, \q\ w 0.016 
(gray lines and squares), (b) Examples of transformation 
of a single peak j = — 1 SDS for non-zero P3 = 0.005 
(dashed line, open circles) and non-zero a — 0.008 (solid 
gray line, gray squares) . Solid black line and filled circles 
mark the p$ = a = case, k = —0.004 for all solutions. 
The other parameters as for Fig. 1. 

SDS), while the dashed lines in Fig. 2(b) mark the solu- 
tions existing for j = +1 (i?-typc SDS). The latter ones 
are mostly unstable, because for j = +1 the discrete 
diffraction is not compensated by the normal GVD. How- 
ever, i?-type SDSs are still of some interest because they 
are involved with the bifurcations joining them to the T- 
type SDSs. An example of an i?-type solution is shown 
in Fig. 4(a). We have also found that the SDS branches 
shown in Fig. 2 can be smoothly continued using the 
Newton method into the region of nonzero higher order 
dispersions and nonzero values of a, which results in the 
appearance of asymmetries in their profiles, sec Fig. 4(b). 
Interestingly, this creates an opportunity to compensate 
asymmetries due to unavoidable p3 with controllable a. 
A more detailed discussion of the existence and stability 
of SDSs goes beyond our present scope. 

Summary: We reported families of spectral-discrete 
solitons existing due to the balance between the effective 
diffraction in frequency space induced by Raman coher- 
ence, and the spectral trapping created by the normal 
GVD. The SDS concept has allowed us to study the prob- 
lem of short pulse generation using methods of nonlinear 
dynamics. We demonstrated that spectrally broad, but 
still localized, solutions associated with trains of short- 
pulses are typically stable and that solutions strongly lo- 
calized in frequency space, corresponding to the weakly 
modulated waves, tend to be dynamically unstable. Our 
approaches could be useful in the analysis of other non- 
linear systems which generate regular frequency combs. 
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